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ABSTRACT 

We introduce the GALEX Arecibo SDSS Survey (GASS), an on-going large pro- 
gram that is gathering high quality Hi-line spectra using the Arecibo radio telescope 
for an unbiased sample of ^1000 galaxies with stellar masses greater than 10 10 M Q 
and redshifts 0.025 < z < 0.05, selected from the SDSS spectroscopic and GALEX 
imaging surveys. The galaxies are observed until detected or until a low gas mass 
■ fraction limit (1.5—5%) is reached. This paper presents the first Data Release, DR1, 

consisting of ~20% of the final GASS sample. We use this data set to explore the 
main scaling relations of Hi gas fraction with galaxy structure and NUV-r colour. A 
large fraction (~60%) of the galaxies in our sample are detected in Hi. Even at stellar 
masses above 10 11 M©, the detected fraction does not fall below ^40%. We find that 
the atomic gas fraction Mhi/M* decreases strongly with stellar mass, stellar surface 
mass density and NUV— r colour, but is only weakly correlated with galaxy bulge-to- 
disk ratio (as measured by the concentration index of the r-band light). We also find 
that the fraction of galaxies with significant (more than a few percent) Hi decreases 
sharply above a characteristic stellar surface mass density of 10 85 M Q kpc -2 . The 
fraction of gas-rich galaxies decreases much more smoothly with stellar mass. One of 
the key goals of the GASS survey is to identify and quantify the incidence of galaxies 
that are transitioning between the blue, star-forming cloud and the red sequence of 
passively-evolving galaxies. Likely transition candidates can be identified as outliers 
from the mean scaling relations between Mhi/M* and other galaxy properties. We 
have fit a plane to the 2-dimensional relation between Hi mass fraction, stellar surface 
mass density, and NUV— r colour. Interesting outliers from this plane include gas-rich 
red sequence galaxies that may be in the process of regrowing their disks, as well as 
blue, but gas-poor spirals. 
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1 INTRODUCTION 

There are good reasons to investigate how and why the 
cold gas content of a massive galaxy varies with stellar 
mass and other physical properties relating to its growth 
history. While the distinction between red, old ellipticals 
and blue, star-forming spirals has been known for a long 
time, recen t work based on the Sloan Digital Sky Sur- 
vey (SPSS; [York et all |2000T ) has shown that galaxies ap- 
pear to divide into two distinct "families " at a stellar mass 
M*~3 xlO 10 M m (IStrateva et al] l200ll ; iKauffmann et al.1 



l2003l ; lBaldrv" et al ]|2004l ). Lower mass galaxies typically have 
young stellar populations, low surface mass densities and 
the low concentrations characteristic of disks. On the other 
hand, galaxies with old stellar populations, high surface 
mass densities and the high concentrations typical of bulges 
tend to have higher mass. New theoretical work has led to a 
diverse set of possible mechanisms to explain this character- 
istic mass scale where galaxies transition from young to old 
llKeres et al.1 120051 : IPekel fc Birnboiml 120061 : iHopkins etaD 
2008), with nearly all o perating via quench ing or regulation 
of the gas supply (e.q.. lMartin et aLll2007l ). Observations of 
the cold Hi gas component - the source of the material that 
will eventually form stars - in galaxies across the transition 
mass, will provide an important new test of those models. 

Initial clues can come from the study of the Hi scaling 
relations of massive galaxi es. Early seminal work (e.g., 
Havnes fc Giovanellil Il984l ; iKnapp. Turner, fc Cunniffe 



19851 ; IRoberts et al.l Il99ll ; see iRoberts fc Havnes! Il994 
for a review) took the first step towards answering 
the question of how the Hi properties of galaxies vary 
as a function of morphological type, environment and 
other physical properties. Although massive "transition" 
galaxies are f ound in these and more recent samples 
llPaturel et al.1 12003!; ISpringob et al.1 120051 . hereafter S05; 
iBothwell. Kennicutt fc Led 1200^ 7! the selection criteria 
make it difficult to identify these in a robust way. One needs 
to quantify the distribution of gas fraction as a function of 
M*, colour, star formation rate (SFR) and other galaxy 
properties in order to understand which galaxies are more 
gas rich or gas poor than the mean. 

It has become common practic e to use "photo- 
metric gas-fractio ns" (IBell et al.l 120031 ; iKannappanf |2004 
IZhang et all 120091 '). which exploit the well-known connec- 
tion b etween SFR and gas content ( Schmidt 1959; Kcn nicuttl 
Il998l ). as a substitute for real gas measurements. However, 
these derived average relations cannot be used to study if 
and how the gas content relates to other properties and phys- 
ical conditions in the galaxies. 

Hi studies of transition objects require large and uni- 
form samples spanning a wide range in gas fraction, stel- 
lar mass and other galaxy properties {e.g., structural pa- 
rameters and star formation). Although blind surveys offer 
the required uniformity, Hi studies of transition galaxies are 
currently not possible because the depths reached by exist- 
ing wide-area blind Hi surveys are very shallow compared 
to s urveys such as the SPSS. The Hi Pa rkes All-Sky Sur- 
vey (|Barnes et al.ll200ll ; lMever et al.ll20oi) . covered ~ 30000 
deg 2 and produced a final catalog of around 5000 Hi detec- 
tions with a median redshift of 2800 km s _ . This should be 
contrasted with the main SDSS spectroscopic survey, which 
covers around 7000 deg 2 and contains more than half a mil- 



lion galaxies with a median redshift of 30,000 km s . The 
recently initiated Arecibo Leg acy Fast ALFA survey (AL- 
FALFA ; [Gi^nenriral][2003) is mapping 7000 deg 2 to con- 
siderably deeper limits. With a median redshift of ~9100 
km s - , ALFALFA for the first time samples the Hi popu- 
lation over a cosmologically fair volume, and is expected to 
detect ~30,000 extragalactic Hl-line sources out to redshifts 
of z ~ 0.06. Even so, the galaxies in the transition regime 
detected with ALFALFA will be predominantly gas-rich. 

In this paper we describe the first results and data re- 
lease from the GALEX Arecibo SPSS Surve£| (GASS), a 
new Hi survey specifically designed to obtain Hi measure- 
ments of ~1000 galaxies in the local universe (0.025 < z < 
0.05) with stellar masses M 4 > 10 10 Mq . As we discuss be- 
low, we observe these massive galaxies down to a low gas 
mass fraction limit (1.5—5%), in order to study the physical 
mechanisms that shape the stellar mass function, regulate 
gas accretion and quench further galaxy growth by conver- 
sion of gas into stars. We expect that GASS will provide a 
rich, homogeneous data set of structural and physical pa- 
rameters (e.g., luminosity, stellar mass, size, surface bright- 
ness, gas-phase and stellar metallicities, AGN content, ve- 
locity dispersion), star formation rates and gas properties. 
Analysis of this unique sample should allow us, for the first 
time, to investigate how the cold gas responds to a variety 
of different physical conditions in the galaxy and obtain new 
insights on the physical processes responsible for the tran- 
sition between blue, star-forming spirals and red, passively- 
evolving ellipticals. 

In a companion paper (Schiminovich et al., in prepara- 
tion; hereafter Paper II), we derive volume- averaged quanti- 
ties for our GASS sample to determine the relative fraction 
of Hi associated with massive galaxies in the local universe, 
and compare with the SFR density to explore how the gas 
consumption timescale varies across the galaxy population. 
In both papers we also discuss how we expect to refine and 
improve our analyses using the full GASS data set. 

In this first paper we describe GASS survey design and 
selection criteria (§ [5] and § [3]). Arecibo observations and 
data processing are discussed in § [4] We provide catalogs of 
SPSS/GALEX parameters and Hl-line spectroscopy mea- 
surements for the 176 galaxies in this first Pata Release 
(PR1) in §[5]and §|6] respectively. Our results are presented 
in § [7] We characterize the properties of the PR1 data set 
in § 17.11 In order to obtain a sample that is unbiased in 
terms of Hi properties, we need to correct for the fact that 
we do not re-observe objects already detected by ALFALFA 
or galaxies found in the Cornell Hi archive (S05). We con- 
struct such a representative sample in § 17.21 and we use it to 
quantify average gas fraction scaling relations as a function 
of other galaxy parameters in § 17.31 and § 17.41 Our findings 
are summarized and further discussed in § [8] 

All the distance-dependent quantities in this work are 
computed assuming £1 — 0.3, A = 0.7 and Ho = 70 
km s~ x Mpc~\ AB magnitudes are used throughout the 
paper. 



http://www.rnpa-garching.mpg.de/GASS/ 
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2 SURVEY DESIGN 

GASS is designed to efficiently measure the Hi content of an 
unbiased sample of ~1000 massive galaxies, for which SDSS 
spectroscopy and GALEX l|Martin et al.|[2005h imaging are 
also available. As described below, the targets are selected 
only by redshift and stellar mass, and observed with the 
Arecibo radio telescope until detected or until a gas frac- 
tion limit of 1.5 — 5% is reached (i.e., a gas fraction limit an 
order of magnitude lower than in objects of similar stellar 
mass detected by ALFALFA at the same redshifts). We de- 
scribe below our survey requirements and sample selection 
methodology. 

Survey footprint. All the GASS targets are located 
within the intersection of the footprints of the SDSS pri- 
mary spectroscopic survey, the projected GALEX Medium 
Imaging Survey (MIS) and ALFALFA. The SDSS primary 
spectroscopic sample targets all galaxies with r < 17.77 with 
high completeness (> 80% for r > 14.5). The GALEX MIS 
reaches limiting NUV magnitude ~23, which allows us to 
probe the full range of colours (and derived SFRs) of the 
normal galaxy population. Existing ALFALFA coverage in- 
creases our survey efficiency by allowing us to remove from 
the GASS target list any objects already detected by AL- 
FALFA — this amounts to an estimated 20% of the galaxies 
meeting our selection criteria. However, this does not cor- 
respond to a 20% gain in observing time, because the ob- 
jects that we skip are those that would be detected with 
the shortest integrations. We also do not re-observe galaxies 
with detections in the Cornell Hi digital archive of targeted 
observations (S05), a homogeneous compilation of Hi param- 
eters for ~9000 optically-selected galaxies, mostly selected 
for Tully-Fisher applications. 

Stellar mass range (10 < LogM*/M Q < 11.5). We 
target a stellar mass range that straddles the "transition 
mass" at ~3 x 10 10 M . 

Redshift range (0.025 < z < 0.05). An Hi sur- 
vey of massive galaxies selected from SDSS is ideally per- 
formed at redshifts above z > 0.025. At magnitudes brighter 
than r ~ 13 (corresponding to a transition mass galaxy 
at z = 0.025), the spectroscopic targeting completeness in 
the SDSS falls below 50%. Additionally, a single pointing 
on galaxies at z < 0.025 may occasionally underestimate 
the Hi flux, if their Hi disks are extended in comparison 
with the Arecibo beam. The upper end of our redshift in- 
terval is set by practical sensitivity limits as well as a desire 
to remain within the velocity range covered by ALFALFA 
(0 < z < 0.06). We further restricted the range to z < 0.05 
in order to avoid the gap in ALFALFA velocity coverage 
caused by radio frequency interference (RFI) at 1350 MHz 
from the the Federal Aviation Administration (FAA) radar 
in San Juan. 

Gas mass fraction/gas mass limit. A crucial goal 
of GASS is to identify galaxies that show signs of recent ac- 
cretion and/or quenching. We wish to obtain accurate gas 
mass measurements for transition galaxies, which have had a 
small, but significant amount of recent star formation (1—5% 
of their total mass). This translates into a requirement that 
we observe the sample to an equivalent gas mass fraction 
(defined as Mhi/M* in this work) limit. Practically, we have 
set a limit of M H i/M* > 0.015 for galaxies with M+> 10 10 ' 5 
Mq, and a constant gas mass limit Mhi= 10 8 ' 7 Mq for 



galaxies with smaller stellar masses. This corresponds to a 
gas fraction limit 0.015 — 0.05 for the whole sample. This 
allows us to detect galaxies with gas fractions significantly 
below those of the Hi-rich ALFALFA detections at the same 
redshifts, and find early-type transition galaxies harboring 
significant reservoirs of gas. We do not try to detect incon- 
sequential amounts of gas (Mhi/M* < 0.01) typical of the 
most gas-poor early-types. 

Based on the Hi mass limit assigned to each galaxy (i.e., 
Mhi= 10 8 ' 7 Mq or 0.015 M*, whichever is larger), we have 
computed the observing time, T max , required to reach that 
value with our observing mode and instrumental setup (see 

§0. 

3 SAMPLE SELECTION 

Since the ALFALFA and GALEX surveys are on-going, 
we h ave defined a GASS parent sa mple, based on SDSS 
DR6 l|Adelman-McCarthv etal]|2008l ) and the maximal AL- 
FALFA footprint, from which the targets for Arecibo ob- 
servations are extracted. The parent sample includes 12006 
galaxies that satisfy our stellar mass and redshift selection 
criteria (see Fig. [I); of these, ~ 10000 have UV photometry 
from either the GALEX All-sky Imaging Survey (AIS; ~100 
s exposure, FUV, im , NUVfcm < 21 m AB ) or the MIS (~1500 
s exposure, FUV iim , NUVjim < 23 m AB ). The final GASS 
sample will include ~1000 galaxies, chosen by randomly se- 
lecting a subset which balances the distribution across stel- 
lar mass and which maximizes existing GALEX exposure 
time. In practice, we have extracted a subset of few hundred 
targets distributed across the whole GASS footprint, which 
includes galaxies that already have MIS or at least AIS data, 
and we have given highest priority to those in the sky re- 
gions already catalogued by ALFALFA. We have also given 
some priority to objects with stellar mass greater than 10 10 ' 5 
Mq, but the final survey sample will have similar numbers 
of galaxies in each stellar mass bin. 

3.1 Overlap with ALFALFA and Cornell Digital 
Hi archive 

As mentioned in §[21 we do not re-observe galaxies with good 
Hi measurements already available from either ALFALFA 
or the S05 digital archive. The overlap between the GASS 
parent sample and the S05 compilation of detected galaxies 
(their Table 3) is 430 objects. ALFALFA has released five 
catalogs to date, three of which are in sky region s with SDSS 
spectroscopic coverage |Giovanelli et al.l 120071 ; iKent et ail 
l200Sl ; IStierwalt et alj|2009l ). Together with still unpublished 
data, ALFALFA has fully catalogued the following two sky 
regions relevant for GASS: (a) 7.5 < «20oo < 16.5 hrs, 
+4° < <5 2 ooo < +16°, and (b) 22 < a 2 ooo < 3 hrs, 
+ 14° < (5 2 ooo < +16°. Notice that the only ALFALFA de- 
tections that we do not re-observe are the ones classified 
as reliable (i.e., "code 1" in their tables); we do target the 
ALFALFA priors, i.e. candidate sources with lower signal- 
to-noise (< 6.5) but optical counterparts with known and 
matching redshift ("code 2"). These sources are usually con- 
firmed by GASS observations and detected with short inte- 
grations (13 of the 176 DR1 galaxies are classified as AL- 
FALFA code 2. They were all detected except one). In the 
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Figure 1. Sky footprint of the GASS survey. Gray points represent galaxies in the parent sample, the super-set of ~12000 objects out 
of which targets for Arecibo observations are extracted. Blue points indicate the subset of galaxies detected by the ALFALFA survey to 
date. 



rest of this paper, "ALFALFA detections" will refer only to 
the galaxies in the first category. Figure [1] shows the sky dis- 
tribution of the 769 sources meeting GASS selection criteria 
that have been detected by the ALFALFA survey to date 
(blue points). These are divided into 658 in the region (a) 
defined above and 81 in (b), corresponding to 18% and 19% 
of the parent sample in the same areas. The S05 archive con- 
tributes 195 galaxies to region (a) and 45 to (b), the large 
majority of which are also ALFALFA detections (but the 
S05 Hi profiles have higher signal-to noise). After removing 
duplicates, the numbers of available Hi detections in regions 
(a) and (b) are 710 and 89, respectively. 



4 ARECIBO OBSERVATIONS AND DATA 
REDUCTION 

GASS observations started in March 2008 and are on-going. 
In the first 1.5 years of the survey until end of August 2009, 
Arecibo allocated 239 hours to this project, with approxi- 
mately 40 hours lost to technical or RFI problems. The ob- 
servations were scheduled in 85 blocks of 1—6.25 hours, with 
about 50% of the total allocation time in blocks of 3 hours 
or less. The full survey will require a total of ~900 hours to 
be completed (not including the increased overhead caused 
by the scheduling in short blocks). Most of the observing 
(~160 hours) has been carried out remotely from the MPA 
in Garching. Since May 2009, remote observations are also 
done from Columbia, JHU, and NYU. 

The Hi observations are done in standard position- 
switching mode: each observation consists of an on/off 
source pair, each typically integrated for 5 minutes, followed 
by the firing of a calibration noise diode. We use the L- 
band wide receiver, which operates in the frequency range 
1120-1730 MHz, with a 1280-1470 MHz filter to limit the 
impact of RFI on our observations. The interim correla- 
tor is used as a backend. The spectra are recorded every 
second with 9-level sampling. Two correlator boards, each 
configured for 12.5 MHz bandwidth, one polarization, and 
2048 channels per spectrum (yielding a velocity resolution 



of 1.4 km s _1 at 1370 MHz before smoothing) are centered 
at or near the frequency corresponding to the SDSS redshift 
of the target. The other two boards are configured for 25 
MHz bandwidth, two polarizations, 1024 channels per spec- 
trum, and are centered at 1365 and 1385 MHz, respectively. 
This setup allows us to monitor the full frequency interval 
of the GASS targets (1353 to 1386 MHz, corresponding to 
z — 0.050 and z = 0.025, respectively) for RFI and other 
problems. All the observations are done during night-time to 
minimize the impact of RFI and solar standing waves on our 
data. The Doppler correction for the motion of the Earth is 
applied during off-line processing. 

A radar blanker is used for ~2/3 of our targets (those 
below 1375 MHz) to avoid RFI caused by harmonics of the 
FAA airport radar in San Juan, which transmits at 1330 
and 1350 MHz. This device is synchronized with the pulsed 
signal of the radar — the data acquisition with the interim 
correlator is effectively interrupted for a time St during and 
after each radar pulse (5t can be chosen by the observer to 
lie between 100 and 750 /js, and is typically set to 400 /is to 
block, in addition to the radar pulse, its delayed reflections 
from aircraft or other obstacles). 

The target lists for the observing runs are prepared 
in advance. The selection is made from the compilation of 
galaxies mentioned in § [3] only a small subset of which are 
visible during each given observing window. The GASS tar- 
gets are randomly chosen from that subset, after excluding 
galaxies with prior Hi detections from ALFALFA or the S05 
archive, or with strong continuum sources within the beam 
that would cause ripples in the baselines. 

We usually acquire up to three 5-minute on-off pairs per 
object per observing session, and accumulate observations 
from multiple sessions. One such pair requires approximately 
13.5 or 11.3 minutes with or without the radar blanker, re- 
spectively, and exclusive of slew-to-source time. For galaxies 
with T max < 4 minutes, or to use end-of-run blocks of or- 
der of 10 minutes, we acquire 4-minute pairs instead, but 
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not shortei[f|. This is a good compromise between obtain- 
ing good quality profiles with small time investment for Hl- 
rich objects and minimizing overheads (which increase when 
an observation is broken into smaller time segments). On- 
source integration times for the sample presented in this 
work ranged between 4 and 90 minutes, with an average of 
13 minutes for the detections and 23 minutes for the non- 
detections (total times are ~2.5 times longer). 

GASS spectra are quickly combined and processed 
during the observations in order to assess data quality, and 
to allow us to stop the integration when the object is de- 
tected, or when the maximum time to reach its limiting gas 
fraction has been reached. A more careful data reduction, 
which includes RFI excision, is performed off-line at a later 
stage. 

The data reduction is performed in the IDL environ- 
ment using our own routines, which are based on the stan- 
dard Arecibo data processing library developed by Phil 
Perillat. More specifically, we adapted the software writ- 
ten to process and m easure the observations described in 
ICatinella et al.l (|2008l ) . which, apart for targeting higher red- 
shift objects, adopted identical observing mode and setup. 
In summary, the data reduction of each polarization and 
on-off pair includes Hanning smoothing, bandpass subtrac- 
tion, RFI excision, and flux calibration. A total spectrum is 
obtained for each of the two orthogonal linear polarizations 
by combining good quality records (those without serious 
RFI or standing waves). Each pair is weighted by a factor 
1/rrns 2 , where rms is the root mean square noise measured 
in the signal-free portion of the spectrum. The two polar- 
izations are separately inspected (they usually agree well. 
If present, polarization mismatches are noted in the Ap- 
pendix), and averaged to produce the final spectrum. 

After boxcar smoothing and baseline subtraction, the 
Hl-line profiles are ready for the measurement of redshift, 
rotational velocity and integrated Hi line flux. Recessional 
and rotational velocities are measured at the 50% peak 
level from linear fits to the edges of the Hi profile. Our 
measurement technique is expla i ned in more detail, e.g., in 
ICatinella. Havnes. fc Giovanellil l|2007l . §2.2). 



5 SDSS AND GALEX DATA 

This section summarizes the quantities derived from optical 
and UV data used in this paper. All the optical parameters 
listed below were obtained from Structured Query Language 
(SQL) queries to the SDSS DR7 database served, unless 
otherwise noted. 

The GALEX UV photometry for our sam ple was com- 
pletely reprocessed by us, as explained in IWang et al.l 
(2009). Briefly, we registered GALEX NUV/FUV and SDSS 
r-band images, and convolved the latter to the (lower reso- 
lution) UV Point Spread Function (PSF) using Image Re- 



2 We considered using 2-minute integrations for galaxies classi- 
fied as "code 2" detections in ALFALFA (see § 13, It . which has an 
effective integration time of 48 seconds. We targeted two of these 
objects during the very first run of the survey - one was detected 
(but the baseline was less than optimal) and one was not. 

3 http://cas.sdss.org/dr7/en/tools/search/sql.asp 



duction and Analysis Facility (IRAF) tasks. The UV PSFs 
are measured from stacked stellar images, obtained by co- 
adding the stars within 1200 pixels from the center of the 
frame. After masking out nearby sources detected in either 
UV or convolved SDSS r-b and images, we used SExtrac- 
tor ()Bertin fc Arnoutsl[l996l ) to calculate magnitudes within 
Kron elliptical apertures, defined on the convolved SDSS im- 
ages. 

The NUV— r colours thus derived are c orrect ed for 
Galactic extinction following IWvder et al] l|2007h . who 
adopted A(X)/E(B - V) = 2.751 for the SDSS r-band and 
A(X)/E(B - V) = 8.2 for GALEX NUV. From these as- 
sumptions, the correction to be applied to NUV— r colours 
is Anuv — A r = 1.9807j4 r , where the extinction A r is ob- 
tained from the SDSS data base (listed in Table Q] below as 
"ext r "). 

Internal dust attenuation corrections are very uncertain 
for galaxies outside the b lue sequence, especia l ly in absence 
of far infrared data ( e.g., Johnson et aDl2007l : IWvder et al.l 
120071 ; ICortese et~aH I2008T ). Moreover, obtaining reliable 
SFRs from NUV photometry is not trivi al in the poorly cal- 
ibrat e d, low specific SFR regime (e.g., ISchiminovich et al.l 
120071 ; ISalim et al. I l2007h . A detailed discussion of these 
issues is beyond the scope of this work, and is presented 
in Paper II. Hence, we do not correct NUV— r colours for 
dust attenuation, nor we derive dust-corrected SFRs in this 
paper. 

Table Q] lists the relevant SDSS and UV quantities for 
the GASS objects published in this work, ordered by increas- 
ing right ascension: 

Cols. (1) and (2): GASS and SDSS identifiers. 
Col. (3): SDSS redshift, zsdss- The typical uncertainty of 
SDSS redshifts for this sample is 0.0002. 
Col. (4): base-10 logarithm of the stellar mass, M*, in so- 
lar units. Stellar masses are derived fr om SDSS photom e- 
try using the m ethodology described in lSalim et al. I l2007l (a 
IChabrieri[2003l initial mass function is assumed). Over our 
required stellar mass range, these values are believed to be 
accurate to better than 30%, significantly smaller than the 
uncertainty on other derived physical parameters such as 
star formation rates. This accuracy in M t is more than suf- 
ficient for this study. 

Col. (5): radius containing 50% of the Petrosian flux in z- 
band, Rso.z, in arcsec. 

Cols. (6) and (7): radii containing 50% and 90% of the Pet- 
rosian flux in r-band, -R50 and Rgo respectively, in arcsec 
(for brevity, we omit the subscript "r" from these quantities 
throughout the paper). 

Col. (8): base-10 logarithm of the stellar mass surface den- 
sity, /i*, in Mq kpc -2 . This quantity is defined as = 
Mi, / (2-kRIq z ), with i?5o, z in kpc units. 
Col. (9): Galactic extinction in r-band, ext r , in magnitudes, 
from SDSS. 

Col. (10): r-band model magnitude from SDSS, r, corrected 
for Galactic extinction. 

Col. (11): NUV— r observed colour from our reprocessed 

photometry, corrected for Galactic extinction. 

Col. (12): exposure time of GALEX NUV image, Tnuv, in 

seconds. 

Col. (13): maximum on-source integration time, T max , re- 
quired to reach the limiting Hi mass fraction, in minutes (see 
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§ [2]) . Given the Hi mass limit of the galaxy (set by its gas 
fraction limit and stellar mass), we computed the required 
integration time to reach this limit at the galaxy's redshift, 
assuming a 5a signal with 300 km s _1 velocity width and 
the instrumental parameters typical of our observations (i.e., 
gain ~10 K Jy -1 and system temperature ~28 K at 1370 
MHz). 



6 Hi SOURCE CATALOGS 

In this section we present the main Hi parameters of the 
99 galaxies detected by GASS to date, and provide upper 
limits for the 77 objects that were not detected. 

Table [2] lists the derived Hi quantities for the detected 
galaxies (ordered by increasing right ascension), namely: 
Cols. (1) and (2): GASS and SDSS identifiers. 
Col. (3): SDSS redshift, zsdss, repeated here from Table[T]to 
facilitate the comparison with the Hi measurement (col. 6). 
Col. (4): on-source integration time of the Arecibo observa- 
tion, Ton, in minutes. This number refers to on scans that 
were actually combined, and does not account for possible 
losses due to RFI excision (usually negligible). 
Col. (5): velocity resolution of the final, smoothed spectrum 
in km s _1 . 

Col. (6): redshift, z, measured from the Hi spectrum. The 
error on the corresponding heliocentric velocity, cz, is half 
the error on the width, tabulated in the following column. 
Col. (7): observed velocity width of the source line profile in 
km s _1 , Wso, measured at the 50% level of each peak. The 
error on the width is the sum in quadrature of the statisti- 
cal and systematic uncertainties in km s _1 . Statistical errors 
depend primarily on the signal-to-noise of the Hi spectrum, 
and are obtained from the rms noise of the linear fits to the 
edges of the Hi profile. Systematic errors depend on the sub- 
jective choice o f the Hi signal boundar ies, and are estimated 
as explained in lGiovanelli et ail (|2007T ). These are negligible 
for most of the galaxies in this sample (only 17 objects have 
systematic errors greater than zero). 

Col. (8): velocity width corrected for instrumental broad- 
ening and cosmological redshift only, Wso c , in km s" 1 . No 
inclination or turbulent motion corrections are applied. 
Col. (9): observed, integrated Hi- line flux density in Jy 
km s _1 , F = J S dv, measured on the smoothed and 
baseline-subtracted spectrum. The reported uncertainty is 
the sum in quadrature of the statistical and systematic er- 
rors (see col. 7). The statistical errors are calculated accord- 
ing to equation 2 of S05. 

Col. (10): rms noise of the observation in mjy, measured on 
the signal- and RFI-free portion of the smoothed spectrum. 
Col. (11): signal -to-noise ratio of the Hi spectrum, S/N, esti- 
mated following ISaintongel (I2007T ) and adapted to the veloc- 
ity resolution of the spectrum. This is the definition of S /N 
adopted by ALFALFA, which accounts for the fact that for 
the same peak flux a broader spectrum has more signal. 
Col. (12): base-10 logarithm of the Hi mass, Mm, in solar 
units, computed via: 



Mm 

M 



2.356 x 10 5 

1 + 2 



dh(z) 



Mpc 



J S dv 
Jy km s~ 



(1) 



where di,(z) is the luminosity distance to the galaxy at 
redshift z as measured from the Hi spectrum. 



Col. (13): base-10 logarithm of the Hi mass fraction, 
Mm/ 

Col. (14): quality flag, Q (l=good, 2=marginal, 5=con- 
fused). An asterisk indicates the presence of a note for the 
source in the Appendix. Code 1 refers to reliable detections, 
with a S/N ratio of order of 6.5 or higher (this is the same 
threshold adopted by ALFALFA). Marginal detections have 
lower S/N, thus more uncertain Hi parameters, but are 
still secure detections, with Hi redshift consistent with the 
SDSS one. The S/N limit is not strict, but depends also 
on Hi profile and baseline quality. As a result, galaxies 
with S/N slightly above the threshold but with uncertain 
profile or bad baseline may be flagged with a code 2, and 
objects with S/N < 6.5 and Hi profile with well-defined 
edges may be classified as code 1. We assigned the quality 
flag 5 to four "confused" galaxies, where most of the Hi 
emission is believed to come from another source within 
the Arecibo beam. For some of the galaxies, the presence of 
small companions within the beam might contaminate (but 
is unlikely to dominate) the Hi signal - this is just noted in 
the Appendix. Finally, we assigned code 3 to GASS 9463, 
which is both marginal and confused. 

Table [3] gives the derived Hi upper limits for the 
non-detections. Columns (1-4) and (5) are the same as 
columns (1-4) and (10) in Table [2] respectively. Column 
(6) lists the upper limit on the Hi mass in solar units, 
Log Mm, Km, computed assuming a 5 a signal with 300 
km s _1 velocity width, if the spectrum was smoothed to 
150 km s" 1 . Column (7) gives the corresponding upper 
limit on the gas fraction, Log Mm.,Um/M+. An asterisk in 
Column (8) indicates the presence of a note for the galaxy 
in the Appendix. 

Figure [2] shows SDSS images and Hi spectra for the 
galaxies with quality flag 1 in Tabled marginal detections 
and galaxies for which confusion is certain are shown sep- 
arately in Figure [3] and non-detections are presented in 
Figure 2] The objects in these figures are ordered by in- 
creasing GASS number (indicated on the top right corner of 
each spectrum). The SDSS images show a 1 arcmin square 
field, i.e., only the central part of the region sampled by 
the Arecibo beam (the half power full width of the beam 
is ~3.5' at the frequencies of our observations). Therefore, 
companions that might be detected in our spectra typically 
are not visible in the postage stamps - examples include 
GASS 40007, a spectacular pair of blue spirals with 1.3' 
separation and 60 km s _1 velocity difference (marked as 
confused in Table [2] even if the Hi spectrum does not ap- 
pear clearly distorted), and a few non-detections (e.g., GASS 
29090, 40686, and 42156). The Hi spectra are always dis- 
played over a 3000 km s _1 velocity interval, which includes 
the full 12.5 MHz bandwidth adopted for our observations. 
The Hl-line profiles are calibrated, smoothed (to a velocity 
resolution between 5 and 21 km s _1 for the detections, as 
listed in Table [2] or to ~15 km s _1 for the non-detections), 
and baseline-subtracted. A red, dotted line indicates the he- 
liocentric velocity corresponding to the optical redshift from 
SDSS. There is a very good agreement between SDSS and 
Hi redshifts, with only small offsets that are usually within 
the typical SDSS measurement uncertainty (0.0002). In Fig- 
ures [2] and [3] the shaded area and two vertical dashes show 
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Figure 2. SDSS postage stamp images (1 arcmin square) and Hl-line profiles of the survey detections, ordered by increasing GASS 
number (indicated in each spectrum). The Hi spectra are calibrated, smoothed and baseline-subtracted. A dotted line and two dashes 
indicate the heliocentric velocity corresponding to the SDSS redshift and the two peaks used for width measurement, respectively. [See 
the electronic edition of the Journal for the complete figure.] 



the part of the profile that was integrated to measure the Hi 
flux and the peaks used for width measurement, respectively. 

The GASS Hi spectral data products will be incorpo- 
rated into the Cornell Hi digital archival, a registered Vir- 



tual Observatory node that already contains the ALFALFA 
data releases and the S05 Hi archive of targeted observa- 
tions. 



http:/ / arecibo. tc. Cornell, edu/hiarchive 
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[See the electronic edition of the Journal for the complete figure.] 
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GASS 9463 is both marginal and confused) detections. 



7 RESULTS 

7.1 GASS DR1 Sample 

We describe here the properties of the galaxies included in 
this first GASS data release, whose Arecibo Hl-line profiles 
and derived parameters were presented in the previous sec- 
tion. 

The sky distribution of the galaxies is shown in Fig- 



ure [5] As can be seen, most of the observing time thus far 
has been allocated in the 110° < «2ooo < 250° region (re- 
ferred to as the "Spring sky", because visible from Arecibo 
during the night-time in that season). The targets are also 
concentrated in the part of sky already covered and cata- 
logued by ALFALFA (rectangles). We did observe several 
galaxies located outside the current ALFALFA catalogued 
footprint, partly because of time allocation constraints and 
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Figure 4. Same as Fig. [2] for non-detections. [See the electronic edition of the Journal for the complete figure] 



partly because, as already mentioned, we gave some priority 



to objects with stellar mass larger than 10 



M 



0- 



The distributions of measured Hi properties for GASS 
detections are presented in Figure [5] (solid histograms) . In 
the top left panel, we show the redshift histogram for the 
full DR1 sample using SDSS measurements (dotted). The 
distribution of velocity widths (not deprojected to edge-on 
view) peaks near 300 km s _1 . As mentioned in the previ- 
ous section, this is the value we adopt to calculate Hi mass 



limits for the non-detections. The measured Hi masses vary 
between 4.6 x 10 s and 3.2 x 10 10 M . Approximately half 
of the galaxies detected by GASS have Hi mass fractions 
smaller than 10%. 

An overview of some of the optical- and UV-derived 
parameters for this data set is found in Figure [7] where 
solid and hatched histograms represent full sample and non- 
detections, respectively (dotted histograms will be discussed 
in the next section). Galaxies that were not detected in Hi 
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Figure 5. Sky distribution of the 176 galaxies in the GASS Data Release 1. The rectangles indicate regions already catalogued by 
ALFALFA to date (ALFALFA coverage of areas without SDSS spectroscopy is not shown). 
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Figure 6. Distributions of redshifts, velocity widths, Hi masses and gas mass fractions for the 99 galaxies with Hi detections. The dotted 
histogram in the top left panel shows the distribution of SDSS redshifts for the full DR1 sample (i.e., including the non-detections). 



are spread throughout the entire GASS stellar mass interval 
(panel a), but are preferentially found at higher M* val- 
ues. The detection rate of the survey as a function of M* 
is shown in panel (b). The fraction of detected galaxies de- 
creases with stellar mass, but does not drop below 35% even 
at stellar masses larger than 10 11 Mq (except for the very 



last bin, where we targeted only two objects). Thus, GASS 
is sensitive enough to detect Hi in a significant fraction 
of massive systems. The bottom panels of Figure [7] show 
SDSS r-band magnitude and observed (i.e., not corrected 
for dust attenuation) NUV— r colour distributions for this 
sample. The well-known separation between blue cloud and 
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Figure 7. Stellar mass (a), SDSS r-band (c) and observed NUV— r colour (d) distributions for the GASS DR1 sample (solid). Hatched 
histograms indicate the corresponding distributions for the non-detections. The detection fraction (i.e., the ratio of detections to total) 
is shown as a function of stellar mass in (b). The dotted histograms in the top panels are the average distributions for the representative 
sample discussed in the text (§ 17.21 1. 



red sequence galaxies, best appreciated w hen UV-to-optical 
colours are used (e.g.. lWvder et aL 2007), is clearly seen m 
panel (d) . Not surprisingly, nearly all the non-detections are 
found in the red sequence, whereas bluer, star-forming ob- 
jects are almost always detected. Particularly interesting are 
the galaxies detected in red sequence. Of the 17 detections 
with NUV— r> 5, half are highly inclined disks, thus their 
colours are likely reddened by dust, but most of the others 
have a featureless, spheroidal appearance in the SDSS im- 
ages. The gas fractions of these detections are all below 6%, 
with two very notable exceptions, GASS 9863 and 3505 (Hi 
mass fractions of 27% and 50%, respectively), both early- 
type galaxies. GASS 3505 in particular is an extraordinary 
system that will be mentioned again in this work — its SDSS 
image and Arecibo detection can be seen in the left column 
of the second row of Figure [5] 

In Figure |8j we examine how the gas mass fraction of 
this sample varies as a function of several other galaxy prop- 
erties. From top to bottom and from left to right, the Hl-to- 
stellar mass ratio is plotted as a function of stellar mass, stel- 
lar mass density, concentration index, and observed NUV— r 
colour. In all panels, GASS detections are represented by 
red circles, and upper limits for non-detections are plotted 
as green, upside-down triangles. For comparison, the Hi-rich 



ALFALFA detection^] in the GASS parent sample are shown 
as small gray circles. Blue circles are randomly-selected, gas- 
rich galaxies from either ALFALFA or the S05 archive of 
pointed observations, added to the DR1 sample in the right 
proportion to quantify average trends in the data (see § 17. 2[) . 
The stars identify the locations of GASS 3505 (red) and 
GASS 7050 (green), which are somewhat extreme examples 
of two opposite types of transition galaxies. The former is 
an early-type galaxy with unusually high gas content that 
might be reaccreting a disk (thus moving from the red to- 
ward the blue sequence), and the latter is a gas-poor, disk 
galaxy (which is likely moving in the opposite direction, from 
the blue to the red sequence. The SDSS image and Hi spec- 
trum of GASS 7050 can be seen in the right column of the 
second row of Figure [4}. This figure illustrates how GASS is 
able to measure the Hi content of massive galaxies down to a 
gas fraction limit which is an order of magnitude lower than 
that achieved by ALFALFA over the same redshift range. 
This allows us to quantify the distribution of atomic gas 
fractions over a much wider dynamic range. The figure also 



5 The Hi masses of galaxies in the ALFALFA and S05 archives 
have been recomputed from the tabulated fluxes using Eq.[T]and 
adopting the cosmological parameters listed in §[T] 
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Figure 8. The Hi mass fraction of the GASS sample is plotted here as a function of stellar mass, stellar mass surface density, concentration 
index, and observed NUV-r colour. Red circles and green triangles represent detections and non-detections, respectively. Blue symbols 
indicate randomly selected, Hi-rich ALFALFA and S05 galaxies added to the DR1 sample in the correct proportion, to compensate for the 
fact that we do not re-observe objects with good detections in either of those archives. The sample including red, green and blue symbols 
is one example of the catalogs described in § I7.2I and contains 193 galaxies. For comparison, we also show the full set of ALFALFA 
galaxies meeting the GASS selection criteria that have been catalogued to date (gray). The stars in all plots highlight the locations of 
two interesting objects, GASS 3505 (red) and GASS 7050 (green), discussed in the text. The dashed line on the top-left panel indicates 
the Hi detection limit of the GASS survey. Dotted lines in each panel are linear fits to the red and blue circles only. 



illustrates the relative strengths of the correlations (or lack 
thereof) between gas content and galaxy structural param- 
eters, and also between gas content and NUV— r, a quantity 
that characterizes the stellar populations of the galaxies in 
our sample. 

We discuss each panel below. In order to quantify 
the scatter in each correlation, we have performed a least- 
squares fit to all the galaxies with Hi detections (i.e., all the 
points plotted as red and blue circles in Fig. [8]). The result 
of the fit is shown as a dotted line in each panel. We only 
report the rms variance in Log Mhi/M* about each relation, 
a, in the text below. 

Stellar mass. — The apparently clean correlation be- 
tween gas mass fraction and M* exhibited by the Hl-rich 
ALFALFA galaxies in this stellar mass regime is mostly due 
to the sensitivity limits of the blind Hi survey. At each stellar 
mass, GASS detects many objects with significantly lower 



gas fractions. The gas fraction of the GASS detections de- 
creases as a function of stellar mass, but the correlation has 
larger scatter (a = 0.386 dex) compared to the one ob- 
tained for the ALFALFA galaxies. As already noted, the 
non-detections span the full M* interval studied by GASS, 
but their fraction of the targeted sample increases with stel- 
lar mass. The stellar mass-dependent detection limit of our 
survey is indicated by a dashed lin^f). 



6 Several data points scatter below the nominal detection thresh- 
old, and some non-detections lie slightly above it. There are two 
main reasons for this. First, the dashed line is the expected limit, 
computed assuming an average value for the telescope gain and a 
5 a signal of fixed velocity width (300 km s" 1 , which is represen- 
tative for the massive galaxies that we target). For the same Hi 
flux, galaxies with narrower profiles (i.e., smaller observed velocity 
widths, which might be intrinsic and/or due to small inclination 
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Stellar mass density. — Interestingly, the gas content 
of massive galaxies seems to correlate better with fi+ than 
with M^. Quantitatively, the rms variance in Log Mm/ M* 
at a fixed value of /x* is a = 0.364 dex, i.e. 6% smaller than 
the variance at a fixed value of M*. Another striking feature 
in this plot is the distribution of the non-detections, which 
all have [i+> 10 8 ' 6 Mq kpc -2 , without a single exception. 
We will come back to this point later. 

Concentration index. — The concentration index is 
defined as R90/R50, where -R90 and R50 are the radii en- 
closing 90% and 50% of the r-band Petro sian flux, respec- 
tively (see Table[TJ . As shown in Figure 1 of lWeinmann et aU 
(2009), there is a tight and well-defined correlation between 
the concentration index and the bulge-to-total ratio de- 
rived from full 2-dim ensional multi-c omponent fits using the 
methods presented in lGadottH (|2009h . It is thus very interest- 
ing that the dependence of Mm/M* ratio on concentration 
is considerably weaker [a = 0.449 dex) than its dependence 
on stellar mass or stellar surface density. Because the Hi 
gas is expected to be located in the disk and not the bulge, 
this lack of correlation might imply that the bulge has little 
influence on the formation or fuelling of the disk. 

NUV— r colour. — Colour is well known to be a rea- 
sonably good predictor of gas content for blue-sequence, 
star-forming galaxies. We show here that the correlation be- 
tween those two quantities continues with increased scatter 
beyond the blue sequence traced by the Hi-rich ALFALFA 
galaxies, and down to the survey gas fra ction limit. This is 
in agr eement with previous work by, e. o.. lCortese fc Hughes! 
( 2009), based on a sample of galaxies located in regions in 
and around nearby clusters. Among the four parameters con- 
sidered here {i.e., M*, /i*, R90/R50 and NUV— r), NUV-r 
colour is the one most tightly correlated with gas fraction 
(a = 0.327 dex). As seen in Figure [TJi, almost all the non- 
detections are found in the red sequence. We also notice that 
GASS 3505 and GASS 7050 stand out as clear outliers in this 
plot, with a gas content well displaced from the average for 
their NUV— r colour. 

In the remainder of the paper, we will quantify the main 
correlations explored in this section and discuss their impli- 
cations. 



7.2 Building a Representative Sample 

Although this first data release amounts to only ~20% of the 
full survey sample, our target selection has not been biased 
by any galaxy property except stellar mass. As described in 
§ O in this initial phase of the survey we have given some 



to the line-of-sight) are easier to detect, thus they might yield 
gas fractions below the dashed line. Non-detections can scatter 
around that line because upper limits are based on the actual 
rms noise measured from the spectra, which might differ slightly 
from the expected value (the telescope gain depends on azimuth, 
zenith angle, and frequency of the observation; the actual rms de- 
pends on baseline quality). Second, we never integrate less than 
4 minutes, even if the maximum time T max computed to reach 
the gas fraction limit is smaller (see §(4j. This affects the higher 
stellar mass galaxies, for which the limit can be reached in as 
little as one minute. As for the detections with small T max val- 
ues, we consider taking an additional on/ off pair when the small 
investment of time guarantees a significantly improved Hi profile. 



priority to galaxies with stellar masses greater than 10 
Mq, but otherwise the selection has been random. We can 
thus use this sample to quantify the correlations between 
Hi mass fraction and other physical parameters discussed 
above, after accounting for the missing Hi-rich objects that 
were not observed because they are in ALFALFA or in the 
S05 archive. 

In order to construct a sample of galaxies that is repre- 
sentative in terms of Hi properties, we clearly need to add 
the correct proportions of ALFALFA and S05 galaxies to the 
GASS data set. We determined such proportions for each 
stellar mass bin as follows. We used the ALFALFA catalog 
footprint available to date to estimate the fractions of GASS 
parent sample galaxies detected by ALFALFA (/a a) and in 
the S05 archive (/sos, not counting those already in /aa). 
One complication is that we observed targets outside that 
footprint, some of which might turn out to be ALFALFA 
detections, and we should not count these twice. Thus, let 
us call JVg and Ng,<aa the total number of GASS objects 
in the given stellar mass bin and those below the ALFALFA 
detection limit, respectively. The numbers of ALFALFA and 
S05 galaxies that should be added back into the sample, for 
each stellar mass bin, are: 



Aaa 



./a 



X /AA - (N G - Ng,<Aa), 



where (Ng — Ag,<aa) effectively count as ALFALFA detec- 
tions (they would be if that survey was already completed), 
and 

AT Ag,<aa , 

JVS05 = Z J X /S05, 

1 - /S05 

where, as noted above, /sos does not include ALFALFA de- 
tections. Repeating this process for each stellar mass bin 
yields a sample of ~200 galaxies. 

In order to better account for statistical fluctuations, 
we generated 100 such catalogs, each time selecting objects 
randomly from the ALFALFA and S05 data sets in the cor- 
rect proportions. More specifically, each catalog is obtained 
by adding N' AA and Ag 05 galaxies to each M* bin of the 
GASS data set, where N' AA and A r g 05 are random Poisson 
deviates of Aaa and Nsos, respectively. The catalogs con- 
tain between 179 and 193 galaxies (187 on average), thus 
the percentage of Hi-rich galaxies added to the GASS sam- 
ple is less than 10%. This is illustrated by the blue symbols 
in Figure [8] for one of these catalogs (red and green symbols 
are GASS detections and non-detections, respectively, and 
blue circles are randomly-selected Hi-rich galaxies that were 
added to each M* bin). 

When we compute the average Hi gas fraction as a func- 
tion of different galaxy properties, we do so based on the 100 
catalogs discussed here. The average data set obtained from 
these catalogs will be referred to henceforth as our represen- 
tative sample. To estimate error bars, we also tried to take 
into account the variance internal to the GASS sample itself 
using standard bootstrapping techniques. We generated 100 
random catalogues in a similar way to that described above, 
but we also draw random indices for GASS galaxies (allow- 
ing for repetitions). The catalogs in this second set include 
between 179 and 200 galaxies each (189 on average). 
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Figure 9. Average trends of Hi mass fraction as a function of stellar mass, stellar mass surface density, concentration index and observed 
NUV-r colour based on the representative sample discussed in 5 17.21 In each panel, large circles indicate average gas fractions (weighted 
as explained in the text). These were computed including the non-detections, whose Hi mass was set to either its upper limit (dark green) 
or to zero (red). Green triangles are weighted medians. The average number of galaxies in each bin, Ni, is indicated above the x axis; 
only bins with TV; > 8 are shown. These results are listed in Table [4] Galaxies in the GASS parent sample detected by ALFALFA are 
plotted as small gray circles. The dashed line in the first panel shows the Hi detection limit of the GASS survey. 



7.3 Gas Fraction Scaling Relations 

In Figure [9j we show how the average Hi mass fraction of 
massive galaxies varies as a function of stellar mass, stel- 
lar mass surface density, concentration index and observed 
NUV— r colour. This quantity is calculated for each of the 
100 catalogs described in the previous section that do not 
include bootstrapping of the GASS galaxies. We average the 
gas fractions in a given bin, properly weighted in order to 
compensate for the uneven stellar mass sampling of the DR1 
data set, using the full GASS parent sample as a reference. 
We bin both the parent sample and the individual catalogs 
by stellar mass (with a 0.2 dex step), and use the ratio be- 
tween the two histograms as a weight. In other words, a 
galaxy in the i-th stellar mass bin and in a given catalog is 
assigned a weight Nps,i/N CSjt ,i, where N P s,i and N ca ,t,i are 
the numbers of objects in the i-th Log bin in the par- 
ent sample and in that catalog, respectively. This is a valid 
procedure because the parent sample is a volume-limited cat- 
alogue of galaxies with M*> 10 10 M©. The final average gas 



fractions are then computed by averaging the weighted mean 
values obtained for each catalog. The results are shown as 
large circles in Figure [9] The difference between green and 
red circles illustrates two different procedures for dealing 
with the non-detections: the Hi mass is set either to the 
upper limit (green) or to zero (red). As can be seen, the an- 
swer is insensitive to the way we treat the galaxies without 
Hi detections, except for the very most massive, dense and 
red galaxies. Note that this is entirely due to the deep limit 
of the GASS survey. The error bar on each circle indicates 
the lcr uncertainty on our estimate of the average gas mass 
fraction, computed from the dispersion of the weighted mean 
values of the bootstrapped catalogs. The average number of 
galaxies that contributed to each bin is indicated in the pan- 
els. Because of the small number statistics in some of these 
bins, our bootstrapping analysis might still underestimate 
the error bars. 
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Figure 10. Gas fraction as a function of stellar mass and stellar mass density for our representative sample. Black symbols and lines are 
weighted averages obtained by setting the non-detections to their upper limits (same as dark green trends in Fig. [9). Left: The sample is 
divided into two bins of stellar mass density (blue: Log fi+< 8.9, red: Log £t*> 8.9). Right: The sample is divided into two bins of stellar 
mass (blue: Log M*< 10.7, red: Log M*> 10.7). Average numbers of galaxies in each bin are indicated in the panels. Only averages 
based on N > 9 galaxies are plotted. 



Weighted mediarQ values of the Mm/M* ratios are plot- 
ted in Figure[5]as triangles. In the highest M*, n+, concentra- 
tion index and NUV— r colour bins, the "median" galaxy is 
a non-detection. The values of weighted average and median 
gas fractions illustrated in this figure are listed in Table|3]for 
reference. Lastly, small gray circles in these panels represent 
galaxies in the GASS parent sample detected by ALFALFA. 
It is clear that the shallower, blind Hi survey is biased to 
significantly higher gas fractions compared to our estimates 
of the global average. 

As the plots in Figure [9] and Table [4] show, the gas 
content of massive galaxies decreases with increasing M*, 
fli, concentration index, and observed NUV— r colour. The 
strongest correlations are with /Lt* and NUV— r. The average 
gas fraction decreases by more than a factor of 30 as in- 
creases from 10 s Mq kpc -2 to a few times 10 9 Mq kpc -2 , i.e. 
the relation is very close to a linear one. A similar large de- 
crease is obtained as a function of NUV— r. In contrast, the 
average Hi fraction decreases by a factor of ~20 over a 1.5 
dex range in stellar mass. The relation with concentration 
index is even shallower, remaining approximately constant 
up to a concentration index of 2.5, and then declining by a 
factor of only ~5 up to the highest values of Rgo/Rso- Fig- 
ure [§] also shows that the difference between the mean and 
median values of Mm/M* is smallest when it is plotted as a 
function of /i* and NUV— r. This is because these two prop- 
erties yield relatively tight correlations without significant 
tails to low values of gas mass fraction (see Fig. [SJ . 

In Figure 1101 we split our sample into two bins in stel- 
lar mass density (left) and in stellar mass (right). Although 
somewhat limited by small number statistics, there is evi- 
dence that a change in Af* has a smaller effect on the re- 



' Given n elements x\...x n with positive weights w\...w n such 
that their sum is 1, the weighted median is denned as the element 
Xk for which: y ' Wi < 1/2 and Wi < 1/2. 



lation between gas fraction and /i* than viceversa. In other 
words, the gas mass fraction of massive galaxies appears to 
be primarily correlated with /i*, and not M*. 

Finally, in Figure 1111 we have plotted the fraction of 
galaxies with Hi gas fractions greater than 0.1 (squares) and 
greater than 0.03 (circles) as a function of stellar mass and 
stellar surface density. As can be seen, the fraction of galax- 
ies with significant (i.e., more than a few percent) gas de- 
creases smoothly as a function of stellar mass. In contrast, 
there appears to be a much more sudden drop in the fraction 
of such galaxies above a characteristic stellar surface density 
of 3 x 10 s Mq kpc -2 . This result can also be seen in the top 
right panel of Figure [8] where we see that all the galaxies 
that were not detected in Hi have stellar surface densities 
greater than 3 x 10 8 M kpc -2 . 

7.4 Predicting the Hi Fraction 

As discussed in the previous sections, the Hi mass fraction is 
most tightly correlated with stellar surface mass density and 
NUV— r colour, and the correlations are close to linear (at 
least for galaxies on the blue sequence). In addition, it can 
be demonstrated that there is also a tight correlation be- 
tween NUV— r colour and stellar surface mass density. This 
implies that a linear combination of the two latter quantities 
might provide us with an excellent predictor for the average 
gas content of a massive galaxy. We have thus fit a plane to 
the 2-dimensional relation between Hi mass fraction, stel- 
lar surface mass densit y, and NUV— r c olour following the 
methodology outlined in lBernardi et all (|2003l ). This is stan- 
dard practice for elliptical galaxies, which have been shown 
to obey a tight plane in the 3-dimensional space of effective 
radius, surface brightness and stellar velocity dispersion. 

We performed this plane-fitting exercise using only de- 
tected galaxies in one of the non-bootstrapped catalogs de- 
scribed in ^ 17.21 f the one illustrated in Figure [SJ. The sample 
used is reproduced in Figure [12] where red and blue circles 
indicate DR1 detections and Hi-rich objects added to the 
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Figure 11. Detection fraction of galaxies as a function of stellar mass (left) and stellar mass surface density (right). Circles and squares 
represent objects with Mm/M* greater that 3% and 10%, respectively. 



sample, respectively. ALFALFA galaxies (small circles) and 
non-detections (upside-down triangles) are not used in the 
fit and are shown for comparison only. The gas fractions 
obtained from our best fit relation are compared with mea- 
sured ones in the figure. The 1:1 relation is indicated by a 
dotted line, and the values of the fit coefficients are given 
in the caption. The rms scatter in Log Mm/ Mi, about this 
relation is 0.315 dex, i.e., we obtain a 4% and 13% reduc- 
tion in the scatter compared to the 1-d relations involving 
NUV— r and /x*, respectively. 

We note that the most gas-rich galaxies in our sample, 
as well as the ALFALFA galaxies, lie systematically above 
the plane, as expected. The most gas-poor galaxies lie sys- 
tematically below the plane. This may result from the fact 
that the non-detections are not used in the fit. With the 
increased sample sizes that will be available to us in future, 
we will be able to stack the non-detections to estimate an 
average gas content, and attempt to use this measurement 
as a way to anchor the prediction at the low gas fraction 
end. 

Perhaps the most interesting use for our best-fit plane 
is as a means to identify interesting objects that deviate 
strongly from the average behavior of the sample. These 
outliers are the best candidates for galaxies that might be 
transitioning between the blue cloud of star-forming spirals 
and the red sequence of passively-evolving galaxies. 

Galaxies which are anomalously gas-rich given their 
colours and densities scatter above the mean relation, while 
those that are gas-poor scatter below. This is clearly demon- 
strated by the Hi-rich ALFALFA galaxies, which are pref- 
erentially found above the line. Also, GASS 3505 (marked 
with a red star on the diagram), a galaxy that has optical 
morphology and colours characteristic of a normal elliptical, 
but a 50% Hi mass fraction, is a clear outlier in this plane. 

Of equal interest are the galaxies with low Hi mass frac- 
tions, but that are still forming stars. These galaxies are 
found near the bottom the plots, but shifted to the right 
(e.g., GASS 7050, a gas-poor disk galaxy that was not de- 
tected in Hi, is indicated by a green star). These may be 
systems where the Hi gas has recently been stripped by tidal 



interactions or by ram-pressure exerted by intergalactic gas, 
or where other feedback processes have expelled the gas. In 
future work, we plan to investigate these different classes of 
transition galaxy in more detail. 



8 SUMMARY AND DISCUSSION 

In this paper we introduce the GALEX Arecibo SDSS Sur- 
vey (GASS), an on-going large program that is gather- 
ing high quality Hi-line spectra for an unbiased sample of 
~1000 massive galaxies using the Arecibo radio telescope. 
The sources are selected from the SDSS spectroscopic and 
GALEX imaging surveys, and have stellar masses M*> 10 10 
Mq and redshifts 0.025 < z < 0.05. They are observed un- 
til detected or until a low gas mass fraction limit (1.5—5% 
depending on stellar mass) is reached. 

Blind Hi surveys such as ALFALFA, which detects 
~20% of the GASS targets, are heavily biased towards 
blue, gas-rich systems at the same redshifts. A number 
of past studies have investigated th e Hi content of ellip- 
tical and early-type galaxies (e.g., M organti et alJ 120061 : 
lOosterloo et al.l [20071 ). but these samples are selected by 
morphology and hence may not be representative of the full 
population of massive galaxies. GASS is the first study to 
specifically target a sample that is homogeneously selected 
by stellar mass, a robust measure that is essential for under- 
standing the observations and connecting them with theory. 

We have presented the first Data Release, DR1, consist- 
ing of ~20% of the final GASS sample. Based on this first 
data installment we have built an unbiased, representative 
sample, which has been used to explore the main scaling 
relations between Hi gas fraction and other parameters re- 
lated to structure and stellar populations of the galaxies in 
this study. Our main findings are discussed below. 

(i) A large fraction (~60%) of the galaxies are detected 
in Hi. Even at stellar masses above 10 11 Mq, the detected 
fraction does not fall below ~40%. Around 9% of galaxies 
more massive than 10 11 Mq have Hi fractions larger than 
0.1 and 34% have Hi fractions larger than 0.03. 
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Figure 12. The best fit "plane" describing the relation between Hi mass fraction, stellar mass surface density and observed NUV-r 
colour is plotted as a dotted line. All the symbols are the same as those in Fig. \E\ The values of the coefficients listed above are the 
following: a = -0.332, b = -0.240, c = 2.856. 



(ii) We have studied the correlation between Mm / M* ra- 
tio and a variety of different galaxy properties. We find that 
the gas fraction of massive galaxies correlates strongly with 
stellar mass, stellar surface mass density and NUV— r colour, 
but correlates only weakly with the concentration index of 
the r-band light. The scatter in the correlations increases 
visibly outside the blue sequence of Hi-rich, star-forming 
spirals. 

(iii) The gas content of massive galaxies decreases for 
increasing values of M* and The latter quantities are 
clearly related, but we presented evidence (although based 
on somewhat limited statistics) that the primary correlation 
is with stellar mass surface density, rather than with stellar 
mass. 

(iv) We also found that the fraction of galaxies with 
significant gas content {i.e., Mhi/M* greater than a few 
percent) decreases strongly above a stellar surface mass 
density of 10 85 M kpc" 2 . This is the threshold stel- 
lar mass density between disk-dominated, late-type galax- 
ies and bulge-dominate d, early-type objects identified by 
iKauffmann et all (|200rJ ) , across which the recent star forma- 
tion histories of local galaxies have been shown to undergo a 
transition. Based on their study of the scatter in colour and 
spectral properties of SDSS galaxies, Kauffmann et al. argue 
that star formation proceeds at the same average rate per 
unit stellar mass below the characteris tic surface density, 
and s huts down above the threshold. ISchiminovich et al.l 
(|2007l ) also noted a striking change in the UV-derived spe- 
cific SFRs of SDSS galaxies at the same characteristic /z+ 
(see their Fig. 20). 

(v) A similar transition in gas prope rties near the charac- 
teristic stellar mass M*~3 xlO 10 Mq (IStrateva et al.ll200ll ; 



IKauffmann et alj 120031 ; iBaldrv et all |2004) is not evident 
from our data. This appears to be in contr adiction to the 
results of iBothwell. Kennicutt fc Led ([2009) , who claim to 
detect such a transition in Mhi/M*, based on a sample of 



much more nearby galaxies. We caution that our results are 
still limited by small statistics, particularly at the low stellar 
mass end. 

One of the key goals of the GASS survey is to identify 
and quantify the incidence of transition objects, which might 
be moving between the blue, star-forming cloud and the red 
sequence of passively-evolving galaxies. Depending on their 
path to or from the red sequence, these objects should show 
signs of recent quenching of star formation or accretion of 
gas, respectively. This task requires us to establish the nor- 
mal gas content of a galaxy of given mass, structural prop- 
erties and star formation rate. The cla ssic concept of Hi de- 
ficiency established for spiral galaxies (IHavnes fc Giovanellil 
11984 ; ISolanes. Giovanelli. fc Havnesl 1 1996f ) is a well-known 
attempt to quantify whether or not a galaxy has been re- 
cently stripped of its Hi gas. Our approach is to fit a plane to 
the 2-dimensional relation between Hi mass fraction, stellar 
surface mass density, and NUV— r colour. We showed pre- 
liminary results that make use of Hi detections only, and 
will include non-detections at a later stage, when we have 
larger samples that will allow us to recover a signal from 
their stacked spectra. 

The GASS survey has already identified a few inter- 
esting examples of transition galaxies. In this work we 
have pointed out two examples, GASS 3505 and GASS 
7050, i.e., a very gas-rich, early- type galaxy and a gas-poor 
disk. Many examples of Hi-deficient galaxies are known 
from the literature. More intriguing are red sequence 
galaxies that might be accreting gas and maybe even 
regrowing a disk (such as perhaps GASS 3505). Two such 
sys tems, NGC 4203 a n d NG C 4262, have been identified 
by ICortese fc Hughes! (|2009l ) from a more local sample 
of galaxies. The authors report convincing evidence that 
the Hi, which is distributed in ring-like structures in both 
NGC objects, might be of external origin. Interestingly, 
also GASS 3505 shows an external ring-like structure 
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in the GALEX images. We have recently obtained VLA 
data for this galaxy, and will report on it in a separate paper. 

The GASS data base of stellar and gas-dynamical mea- 
surements will provide an unprecedented view of the gas 
properties and kinematics of massive galaxies, complement- 
ing the results of on-going blind surveys such as ALFALFA. 
Once completed, GASS will allow us to move beyond the 
mean gas fraction scaling relations studied in this paper, 
and address second order questions, such as how the gas 
fractions depend on metallicity or other quantities at fixed 
stellar mass. GASS will also quantify the frequency with 
which different kinds of transition galaxies occur in the lo- 
cal universe, and how this depends on factors such as local 
environment and AGN content. Insight into the nature of 
massive and transition objects will give us a strong founda- 
tion to further understand the gas properties of high-redshift 
galaxies, and will help guide future directions for Hi surveys 
at existing and planned radio facilities. 
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APPENDIX: NOTES ON INDIVIDUAL 
OBJECTS 

We list here notes for galaxies marked with an asterisk in 
the last column of Tables [2] and [3] The galaxies are ordered 
by increasing GASS number. In what follows, AA1 and 
AA2 are abbreviations for ALFALFA detection codes 1 and 
2, respectively. 

Detections (Tabled]) 

3261 - detected on top of weak RFI across all band; AA1. 
3504 - detected on top of weak RFI across all band; low 
frequency edge uncertain, systematic error. 
3645 - detected on top of weak RFI across all band; AA2. 
3777 - detected on top of RFI across all band; AA2. 
3817 - hints of RFI across band; AA1. 
3962 - disk ~1.5' N has z=0.027, no contamination prob- 
lems; blue spiral ~3' E, z=0.0314 (1377.16 MHz), detected 
in board 4. Low frequency edge uncertain, systematic error. 
3971 - hints of RFI across band in one pair; AA2; low 
frequency edge uncertain, systematic error. 
6375 - group: 2 blue galaxies 2.7' away, GASS 6373 
(= VCC 1016, spiral, similar size) to the S and SDSS 
J122716.82+031812.9 (smaller) to the W, same cz. The two 
spirals do not look distorted and the Hi profile is not clearly 
confused, just slightly offset in z (but notice that GASS 
6373 is offset in the wrong direction) . Contamination likely. 
6506 - GASS 6501, red galaxy ~3' away, 1355.75 MHz not 
detected; RFI spike at 1350 MHz. 
7025 - high frequency peak uncertain. 

7493 - GASS 7457 (non-detection) ~l' away, no con- 
tamination problems (z=0.026; GASS 7457 has z=0.036); 
uncertain profile, systematic error; SDSS emission line z 
(0.02622) in better agreement with Hi. 

7509 - small satellites, some contamination very likely; 
stronger in polarization A; high frequency edge uncertain, 
systematic error. 

9343 - 2 blue galaxies within 1', no cz, possible contami- 
nation; poor fits to edges. 

9463 - merger, plus another blue galaxy ~1' away with 
same z: confusion certain; very uncertain Hi profile, sys- 
tematic error. 

9514 - smaller disk galaxy ~1' away, same z, some con- 
tamination likely; AA2; poor fit to low frequency edge. 
9619 - blue, disrupted little galaxy ~3' away, same z, some 
contamination possible. 

9776 - small blue disk ~15" away, same cz (z=0. 027568, 
1382.30 MHz), some contamination certain. Looks like the 
edges of the Hi profile come from 9776, but the central peak 
is dominated by the companion. 

9814 - blue smudges nearby + small galaxy ~1.5' NW, 
same cz; possible contamination; low frequency peak 
uncertain. 

10019 - RFI spike near 1375.3 MHz, slightly drifting in 
frequency. 

10132 - two small companions: a blue one ~1.3' N, strong 
Ha emission in SDSS spectrum, z=0.0438 (1360.80 MHz), 
and a galaxy ~1.5' E, no Ha in SDSS spectrum, z=0.0444 
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(GASS 10132 has z=0.0443). The strong peak in the Hi 
spectrum is centered on the blue companion. Blend. 
11016 - galaxy pair: 11016 is a red gal, no emission lines in 
SDSS spectrum, centered at 1369.07 MHz, the companion 
is a blue galaxy 1.2' SE, strong emission lines in SDSS 
spectrum, centered at 1368.93 MHz. Confusion certain. 
11223 - small early type ~40" away, no cz, small contami- 
nation possible; poor fit to high frequency edge; AA2. 
11386 - small peak outside Hi profile (on the high fre- 
quency side) is mostly in polarization B. 
11956 companion also detected (SDSS 

J000814.67+150752.9, z=0.0372, 1369.46 MHz, blue 
disk 2.1' SW), no overlap. 

11989 - small blue galaxy 2' away, no cz, but confusion 

unlikely. High frequency edge uncertain, systematic error; 

AA4 (i.e., classified by ALFALFA as tentative detection 

with optical counterpart of unknown redshift). 

12371 - 3 blue galaxies within 3', no contamination 

problems (largest one has z=0.008, other two have z=0.027 

and are detected in board 4 (1383 MHz). 

12983 - low frequency edge uncertain, systematic error. 

13227 - stronger in polarization B; uncertain profile, poor 

fit to high frequency edge, systematic error. 

14831 - small galaxies nearby (within ~1'), no cz; some 

contamination possible; high frequency edge uncertain, 

systematic error. 

15181 - AA1. 

17640 - AA2. 

18279 - confused or blend: most of the signal comes from 
blue galaxy ~1.5' away, z=0.0492 (1353.80 MHz); 2 other 
galaxies within 2' have very different z (0.18 and 0.07). 
18421 - stronger in polarization B. 
18581 - low frequency peak uncertain. 
20133 - galaxy pair: GASS 20165 ~l' away, z=0.0498 
(1353.04 MHz; GASS 20133 has z=0.0489, 1354.19 MHz); 
however 20165 is a red, early type gal, not detected in 10m 
on-source, so significant contamination is unlikely. 
20144 - uncertain profile edges, systematic error. 
20286 - RFI at 1376 MHz. 

23315 - blue disk ~1' away has z=0.054, no contamination; 
low frequency edge uncertain, systematic error. 
23445 - 2 galaxies at ~2.7': one at N has z=0.08, one at 
W is small, has ~same cz (z=0. 047065, 1356.56 MHz) and 
is not detected in spectrum (a bit offset so should be visible 
if there); better in polarization B; uncertain profile. 
24740 - uncertain profile, systematic error. 
26822 - AA2. 

28526 - RFI spike near 1354 MHz (within profile), 2 
channels replaced by interpolation. 

29487 - small companion 2.4' S, z=0.0424 (1362.63 MHz, 
0.5 MHz offset), contamination possible but unlikely; 
strong, narrow RFI near 1359 MHz. 

29505 - much stronger in polarization B; no trace of RFI, 
well centered on SDSS z. 
29842 - AA2. 

30175 - blue companion 1.6' away, slightly lower z 

(z=0.0223, 1389.42 MHz) also detected. Two other galaxies 

~3' away have z=0.08. 

30401 - AA2. 

31156 - AA2. 

38462 - poor fits to edges. 

38703 - stronger in polarization B. 



38717 - narrow RFI at 1360 and 1370 MHz (no radar 
blanker) . 

38758 - small, blue companion ~0.5' away, no cz; galaxy 
clearly distorted, contamination very likely. 
38964 - small satellite to the S? blue, edge-on galaxy > 3' 
away, z=0. 033468 (1.6 MHz away), not detected but might 
be responsible for asymmetry (raising low frequency peak). 
39567 - GASS 39600 ~2.5' away to the N, no contamina- 
tion problems (z=0.044, GASS 39567 has z=0.031); much 
stronger in polarization B. 

39595 - AA2. Low frequency edge uncertain, systematic 
error. 

40007 - galaxy pair (separation 1.3', difference of reces- 
sional velocities is 60 km s _1 ); 2 blue spirals, same cz, ~3' 
and 4.5' away (group); Hi spectrum does not look confused 
but contamination certain (note z off with respect to SDSS 
redshift of both galaxies). 

40024 - group: 3 small galaxies around (l'-3), 2 with 
^sarne cz, one without cz; some contamination likely; high 
frequency edge and peak uncertain, systematic error, poor 
fit. 

40393 - uncertain peaks. 

40494 - RFI spike near 1357 MHz (within profile), 3 
channels replaced by interpolation. 

40781 - strong RFI spike at 1356 MHz (within profile), 4 
channels replaced by interpolation. 

41969 - AA1. 

41970 - small disk galaxy 10" to the W, contamination 
possible (perhaps higher z, though). 

42015 - AA2. 

42167 - 3 disks within 2'-4' have different z (largest one 
at ~4', 1373.17 MHz, not detected), no contamination 
problems. 

47221 - low frequency edge uncertain, systematic error 

(little side peak is in polarization B only). 

47405 - high frequency edge uncertain, systematic error. 

Non-detections (Table [3J 

7286 - small companion ~1.5' away, same z. 

7457 - see GASS 7493 (detection). 

9507 - RFI spike at 1352 MHz. 

9891 - blue disk ~2.5' away, z=0.0243 (1386.71 MHz) not 
detected. 

10150 - galaxy ~1' away has z=0.092, no contamination 
problems. 

10358 - group (3 galaxies ~same cz within 3'). 
10367 - group (3 galaxies ~same cz within 1.5'-4'). 
10404 - small blue galaxy 20" E, no cz, responsible for 
small spike at 1369.3 MHz?. 

12455 - group: 2 other galaxies at 2.5' W and 3.5' E, 
~same cz; marginally detected the one at 2.5' (small 
edge-on disk, z=0.0484, 1354.83 MHz)?. 
12458 - RFI spike at 1370.2 MHz. 

12460 - RFI spike at 1350.8 MHz. Detected blue compan- 
ion (irregular galaxy ~1' away, no cz) at 1352.2 MHz? Not 
seen RFI there, but much stronger in polarization B. 
13156 - AA2. 

13549 - group: several small gals, ~same z, within 3'; also 
a larger early type, no cz. 

18422 - small galaxy ~1.3' W, no cz (+ other, smaller 
ones) . 

21023 - RFI at 1358-59 MHz, slightly drifting in frequency. 
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25154 - galaxy pair, companion is a spiral ~1' away, optical 
em lines, z=0.0364 (1370.52 MHz), also not detected. 
25214 - group (4 galaxies ~same cz within 3'). 
25575 - galaxy pair (galaxy ~same cz 20" W). 
26958 - disturbed morphology, merger? marginally de- 
tected companion (blue disk ~0.5' away, no cz)?. 
29090 - detected blue edge-on disk ~3' S; 3 other disks 
~3'-4' away. 

29420 - group: two small galaxies ~2' away, same cz, + 
others, no cz; looks like a galaxy is detected in the off. 
29699 - group: 2 galaxies ~same z within 2', and another 
one at z=0.044. 

30479 - disk ~1.5' away has z=0.08, no contamination 
problems. 

38472 - blue galaxy ~1' SW has z=0.202, no contamina- 
tion problems. 

39448 - blue galaxy ~1' away, same cz (z=0.0341; GASS 
39448 has z=0.0339), contamination certain. 
39465 - blue irregular galaxy ~1' NW detected in board 
3; galaxy ~1' SW has z=0.076; 2 galaxies ~0.5' SW, no cz, 
marginally detected?. 

39606 - ~3' away from GASS 39607, same cz, also not 
detected. 

39607 - see GASS 39606. 

40257 - group: 3 small galaxies within 2.5' E and a large 
disk ~2' NW, ~same cz. 

40317 - group: 3 galaxies within 2.5', ~same cz; marginally 
detected blue irregular, 2' SW, z=0.0404 (1365.25 MHz)?. 
40686 - group: 4 galaxies within 2', ~same cz, one detected 
(SDSS J131527.88+095243.7, a bluish, edge-on disk -2' 
away, z=0.0505, 1352.12 MHz). 

40790 - small companion ~1' away, z=0.0499 (1352.90 
MHz), also not detected. 

41974 - group in background (3 galaxies at 2'-3' with 
z=0. 045-0. 047), no contamination problems. 
42020 - blue irregular galaxy nearly attached, no cz; neg. 
spike near 1373 MHz is not RFI and is in both pols... galaxy 
in off?. 

42156 - detected SDSS J151719.24+072828.3, a blue disk 
-1.5' SW, z=0.0367 (1370.16 MHz). Possible contribution 
from GASS 42164 at the same cz (1370.03 MHz), a larger, 
face-on, red spiral ~2' away from GASS 42156. Also, small 
disk, 2' E, z=0.077. No confusion problems. 
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Table 1. SDSS and UV Parameters of the GASS Galaxies. 
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Table 2. Hi Properties of GASS Detections. 
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Table 3. GASS Non-detections. 









Ton 




Log M H i Hm 










GASS 




2 SDSS 


(min) 


v II1J y 1 


(M 


3) 


Log M H i , 




Not 


V 1 ) 








(W 

K a ) 


(6) 




f71 




(°7 


12025 


J001934. 54+161215.0 


0.0366 


14 


0.43 


< 9 


10 


< 


-1 


74 




4239 


J015816. 23+141747. 9 


0.0261 


5 


0.58 


< 8 


92 


< 


-1 


88 




14943 


J101600. 20+061505. 2 


0.0458 


4 


0.69 


< 9 


50 


< 


-1 


83 




18900 


J102001. 61+083053. 6 


0.0453 


23 


0.28 


< 9 


09 


< 


-1 


83 




26602 


J103347.41+ 124358.1 


0.0325 


15 


0.40 


< 8 


95 


< 


-1 


80 




23408 


J105322. 36+111050. 4 


0.0430 


10 


0.45 


< 9 


26 


< 


-1 


78 




23450 


J105648.58+120535.7 


0.0476 


48 


0.18 


< 8 


96 


< 


-1 


85 




17659 


J105807.59+091633.9 


0.0344 


4 


0.74 


< 9 


27 


< 


-1 


90 




29699 


J110818. 34+131327.5 


0.0340 


5 


0.66 


< 9 


21 


< 


-1 


80 


* 


12455 


.7112017.79+041913. 3 


0.0491 


4 


0.72 


< 9 


58 


< 


-1 


87 


* 


12460 


J112048. 30+035021.0 


0.0494 


25 


0.29 


< 9 


19 


< 


-1 


79 


* 


12458 


.1112118.26+033953.0 


0.0394 


8 


0.47 


< 9 


20 


< 


-1 


84 


* 


24094 


J122030. 18+112027.3 


0.0431 


10 


0.46 


< 9 


26 


< 


-1 


82 




18422 


J122123. 25+095053.0 


0.0467 


90 


0.17 


< 8 


91 


< 


-1 


65 


* 


12966 


.1123632.24+061010.5 


0.0395 


4 


0.72 


< 9 


38 


< 


-1 


80 




30479 


.1123708.06+142426.9 


0.0308 


70 


0.17 


< 8 


54 


< 


-1 


75 


* 


30471 


J123753. 20+141652. 7 


0.0263 


4 


0.68 


< 9 


00 


< 


-1 


83 




13037 


J124314. 97+040502.0 


0.0485 


18 


0.35 


< 9 


26 


< 


-1 


77 




6565 


J124938. 19+024520.2 


0.0476 


38 


0.28 


< 9 


13 


< 


-1 


76 




6583 


J125055. 79+031149. 3 


0.0483 


8 


0.59 


< 9 


48 


< 


-1 


72 




40500 


J125911. 09+103006.0 


0.0460 


29 


0.28 


< 9 


12 


< 


-1 


76 




30508 


J125926. 22+142030.0 


0.0471 


17 


0.39 


< 9 


28 


< 


-1 


74 




13156 


.1130446.99+035417.8 


0.0408 


9 


0.59 


< 9 


32 


< 


-1 


70 


* 


25154 


J130457.41+120444.6 


0.0358 


4 


0.75 


< 9 


32 


< 


-1 


81 


* 


40570 


J131104.94+084828.3 


0.0325 


4 


0.73 


< 9 


22 


< 


-1 


92 




25214 


J131232. 81 + 114344. 2 


0.0311 


4 


0.68 


< 9 


14 


< 


-2 


02 


* 


40686 


J131529. 82+095100. 6 


0.0496 


10 


0.39 


< 9 


32 


< 


-1 


84 


* 


40790 


J131944.64+102145.9 


0.0489 


9 


0.44 


< 9 


36 


< 


-1 


83 


* 


25347 


J133019. 15+113042. 5 


0.0378 


5 


0.59 


< 9 


26 


< 


-1 


83 




13549 


J134525. 31+034823. 8 


0.0325 


10 


0.49 


< 9 


04 


< 


-1 


78 


* 


7050 


J134909. 69+024511. 5 


0.0489 


23 


0.34 


< 9 


25 


< 


-1 


74 




40317 


J135533. 72+144552. 7 


0.0408 


8 


0.52 


< 9 


27 


< 


-1 


80 


* 


40257 


J135842. 23+132722. 9 


0.0393 




0.63 


< 9 


32 


< 


-1 


81 


* 


40247 


J135942. 61 + 124412. 5 


0.0392 


4 


0.70 


< 9 


37 


< 


-1 


98 




25575 


J140606. 72+123013.6 


0.0379 


5 


0.60 


< 9 


27 


< 


-1 


95 


* 


38529 


J140946. 81 + 113505. 4 


0.0382 


4 


0.63 


< 9 


29 


< 


-2 


03 




7286 


J141432. 05+031124. 9 


0.0269 


10 


0.51 


< 8 


89 


< 


-1 


84 


* 


38472 


J141608. 76+103543.8 


0.0264 


43 


0.22 


< 8 


52 


< 


-1 


67 


* 


38591 


J141740. 51+103459.9 


0.0271 


48 


0.19 


< 8 


46 


< 


-1 


82 




41323 


J141822. 46+080551.0 


0.0440 


13 


0.35 


< 9 


16 


< 


-1 


86 




30811 


J141845. 69+055004. 7 


0.0489 


5 


0.65 


< 9 


53 


< 


-1 


75 




9507 


J142032. 84+050638. 3 


0.0496 


75 


0.18 


< 8 


98 


< 


-1 


75 


* 


7457 


J142713. 78+025048. 6 


0.0357 


15 


0.40 


< 9 


03 


< 


-1 


78 


* 


41482 


J143152. 89+071915.1 


0.0273 


4 


0.74 


< 9 


07 


< 


-2 


00 




29090 


J143810. 20+092009. 7 


0.0303 


4 


0.62 


< 9 


08 


< 


-1 


93 


* 


7581 


.1143915. 52+024340. 9 


0.0279 


4 


0.70 


< 9 


06 


< 


-1 


98 




9748 


J143917.94+032206.0 


0.0279 


s 


0.68 


< 9 


05 


< 


-1 


78 




9917 


.1144025.99+033556.0 


0.0281 


10 


0.46 


< 8 


88 


< 


-1 


78 




9704 


J144059. 30+030813. 5 


0.0265 


4 


0.81 


< 9 


08 


< 


-1 


87 




9891 


J144225. 70+031354. 9 


0.0258 


4 


0.74 


< 9 


02 


< 


-1 


95 


* 


29420 


J144858. 71 + 122924.4 


0.0474 




0.66 


< 9 


51 


< 


-1 


75 


* 


10031 


J145106. 43+045032. 6 


0.0273 


8 


0.55 


< 8 


94 


< 


-1 


89 




10040 


J145235. 22+043245.1 


0.0286 


4 


0.66 


< 9 


06 


< 


-1 


92 




10150 


J150026. 75+041044.4 


0.0328 


34 


0.24 


< 8 


73 


< 


-1 


80 


* 


10292 


J151322. 09+040701. 2 


0.0426 


90 


0.15 


< 8 


77 


< 


-1 


67 




42020 


J151516.46+063918.5 


0.0352 


53 


0.19 


< 8 


71 


< 


-1 


80 


* 


42017 


J151524.84+062654.0 


0.0452 


90 


0.15 


< 8 


84 


< 


-1 


71 




42140 


J151531. 97+072829.0 


0.0457 


18 


0.32 


< 9 


16 


< 


-1 


82 




10367 


J151553. 85+030301.1 


0.0379 


5 


0.67 


< 9 


32 


< 


-1 


76 


* 


10358 


.7151711. 15+032105. 7 


0.0370 


4 


0.78 


< 9 


36 


< 


-1 


73 


* 


10404 


J151722. 96+041248. 9 


0.0361 


5 


0.70 


< 9 


29 


< 


-1 


73 


* 


42156 


.7151724.33+072921. 9 


0.0336 


40 


0.20 


< 8 


67 


< 


-1 


85 


* 


41974 


J151758. 30+064445. 4 


0.0349 


28 


0.24 


< 8 


79 


< 


-1 


85 


* 


42175 


J151832.42+070720.7 


0.0456 


30 


0.28 


<9 


10 


< 


-1 


78 




10447 


J151840. 93+042505. 3 


0.0471 


87 


0.16 


< 8 


88 


< 


-1 


80 




39469 


J151903. 38+080819. 4 


0.0338 


4 


0.84 


< 9 


31 


< 


-1 


81 




39467 


J151953. 53+080557.2 


0.0334 




0.66 


<9 


20 


< 


-1 


78 




39465 


.7152028.70+081706.6 


0.0372 


70 


0.18 


< 8 


73 


< 


-1 


77 


* 


39448 


J152037.22+080305.7 


0.0338 


89 


0.16 


< 8 


59 


< 


-1 


49 


* 


39605 


J152559. 84+094724. 5 


0.0339 


89 


0.15 


<8 


56 


< 


-1 


57 




39607 


J152706. 28+094746. 3 


0.0438 


20 


0.32 


<9 


13 


< 


-1 


77 


* 


39606 


J152716. 73+094603. 8 


0.0437 


69 


0.17 


<8 


84 


< 


-1 


81 


* 


26958 


J154654.33+055328.3 


0.0419 


4 


0.81 


< 9 


49 


< 


-1 


77 


* 


21023 


J155636. 91+272911. 9 


0.0415 


89 


0.16 


< 8 


78 


< 


-1 


54 


* 


10817 


J220120. 93+121148.1 


0.0291 


15 


0.38 


<8 


83 


< 


-1 


78 




10872 


J221321. 50+132611. 3 


0.0281 


24 


0.33 


< 8 


74 


< 


-1 


74 




10884 


J221430. 63+130444. 9 


0.0257 


18 


0.31 


<8 


63 


< 


-1 


84 





© 0000 RAS, MNRAS 000, 000-000 



The GALEX Arecibo SDSS Survey. I. 25 



Table 4. Weighted Average and Median Gas Fractions Plotted in Figure 9 







<M m /M*} 


(M m /M„) 






X 


(*> 


( average) 12 


(average) 6 


(median) c 


{N) d 


ijOg 


10.16 




U.O'iOXU.UUU 


0.206 


24 




10.49 


0.124±0.026 


0.117±0.027 


0.049 


37 




10.77 


0.078±0.023 


0.072±0.024 


0.030 


51 




11.05 


0.050±0.011 


0.042±0.011 


0.017 


60 




11.29 


0.024±0.011 


0.014±0.013 


0.016 


15 


Log /I* 


o.Zo 


n /i o7_|_n ill 


U.^to 1 ZCU.lll 


n 31 r 


i n 

1U 




8.58 


0.236±0.045 


0.235±0.045 


0.149 


20 




8.88 


0.097±0.020 


0.092±0.020 


0.050 


65 




9.15 


0.036±0.007 


0.026±0.007 


0.017 


73 




9.36 


0.020±0.006 


0.011±0.008 


0.017 


12 


-K90 / ttso 


i.oy 


fl 9^^-1-0 H77 


n OQ«_t_n H77 


U. lyo 


1 9 
1Z 




Z.o4 
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n on/i_l_n n^v 


U.U f z 


ZD 




2.70 


0.154±0.024 


0.149±0.024 


0.070 
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0.059±0.011 


0.049±0.012 


0.018 
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0.040±0.015 


0.030±0.017 
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NUV-r 


2.89 


0.280±0.036 


0.280±0.036 


0.208 
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3.60 


0.116±0.016 


0.115±0.016 


0.077 
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4.36 


0.093±0.021 


0.090±0.022 


0.048 


29 




5.33 


0.041±0.008 


0.028±0.008 


0.016 


45 




5.86 


0.019±0.002 


0.006±0.003 


0.016 


41 



Notes. — "Weighted, average gas fraction; Hi mass of non-detections set to upper limit. ^Weighted, average gas fraction; Hi mass of 
non-detections set to zero. c Weighted, median gas fraction; Hi mass of non-detections set to upper limit. d Average number of galaxies 
in bin. 
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